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摘要 : 在 弹簧 作为 航天 任务 中 有 效 载 荷 的 推进 装置 的 优化 设计 中 ,为 保证 载荷 的 释放 速度 以 及 弹 往 质 量 最 轻 , 结 合 几何 约 
束 、 可 靠 性 与 稳定 性 的 要 求 ,建立 非 线 性 约束 优化 模型 ,将 Matlab 求 得 的 理论 解 带 入 商业 软件 Adams 进行 仿真 运算 ,得 到 理 
论 速度 与 仿真 速度 的 拟 合 关系 。 利 用 拟 合 关系 和 用 Matlab 优化 工具 箱 即 可 求 得 满足 仿真 分 离 速度 的 弹簧 设计 参数 。 利 用 
上 述 方法 对 分 离 推力 弹簧 进行 设计 与 计算 , 并 利用 虚拟 样机 技术 进行 验证 ,减少 了 弹簧 设计 重复 试 算 的 过 程 ,弹簧 质量 缩减 
率 在 30% 以 上 ,满足 释放 速度 要 求 , 得 到 满意 的 优化 设计 结果 。 
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Design and Calculation of Compression Spring of Releasing Payload 
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ABSTRACT: Springs can be used to provide release rate for payload in space missions. JIn order to guarantee the 
precise releasing velocity of payload and the minimum mass of springs, a nonlinear constrained optimization model 
was built with geometric constraints, reliability and stability. The theory solution of MATLAB was brought into Adams 
for simulation, which can work out the fitting relationship between theory and simulation. The spring design parame— 


ters were obtained by using the fitting relationship and calculation of MATLAB. The result was verified by Virtual 


Prototype Technology and perfect performance of design parameters was achieved. 
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1 引言 


随 着 科学 技术 的 不 断 发 展 ,人 类 空间 科学 探测 的 研究 领 
域 大 为 拓展 ,空间 物理 探测 已 到 了 多 点 、 多 尺度 的 时 代 , 这 就 
需要 多 颗 小 卫星 或 多 个 有 效 载荷 进行 协同 观测 。 载 荷 分 离 
释放 技术 是 影响 有 效 载荷 初始 状态 和 编队 飞行 的 重要 因素 。 
航天 飞行 器 上 有 效 载荷 的 分 离 机 构 一 般 由 火 工 品 或 形状 记 
忆 合 金 提供 释放 速度 。 但 火 工 品 存 在 安全 性 较 差 ,可 能 会 对 
有 效 载荷 中 的 科学 设备 造成 污染 ,工作 时 产生 较 大 的 冲击 ， 
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只 能 一 次 性 使 用 ,不 能 重复 测试 性 能 等 缺点 ; 而 形状 记忆 合 
金具 有 成 本 较 高 ,作用 时 间 长 ， 需 要 供电 加 热 消耗 大 量 能 源 
等 缺点 ” ,因此 利用 机 械 弹 簧 所 储存 的 弹性 势能 为 低速 有 效 
载荷 提供 初始 速度 就 具有 明显 的 实用 性 和 可 靠 性 。 

弹 咎 作为 一 种 机 械 结构 中 的 常用 部 件 ,在 各 行 各 业 中 广 
泛 利用 。 在 车 辆 工程 机 械 设计 等 研究 领域 中 ,主要 考虑 弹 
簧 的 应 力 、 应 变 、 基 频 和 工作 载荷 等 结构 性 能 参数 ,研究 
方法 主要 有 : 优化 设计 理论 ”、 可 靠 性 优化 理论 "”、 多 目 
标 稳健 优化 设计 " 、 多 目标 模糊 优化 设计 " 、 基 于 区 间 算 
法 的 优化 设计 、 基 于 Matlab 优化 算法 的 优化 设计 中 、 基 
于 有 限 元 法 的 优化 设计 "“-“” , 上述 研 究 均 未 考虑 弹 段 的 动 
力学 性 能 及 运动 过 程 ,而 在 有 效 载荷 ( 以 下 简称 “载荷 ”) 分 
离 释放 技术 中 ,分 离 推力 弹簧 ( 以 下 简称 “弹簧 ") 所 能 提供 


的 分 离 速 度 是 首要 的 性 能 参数 。 对 于 单 载荷 释放 任务 , 文献 
[6] 基于 多 岛 遗 传 算法 对 卫星 的 释放 装置 进行 优化 设计 ,并 
对 其 进行 运动 仿真 与 实验 。 对 于 多 载荷 协同 编队 任务 ,为 保 
证 释放 后 载荷 能 顺利 进行 姿态 控制 和 轨道 控制 , 弹 得 所 提供 
的 释放 速度 必须 为 精确 值 , 且 由 于 安装 空间 有 限 , 弹 签 的 几 
何 尺寸 必须 满足 空间 限制 。 然 而 文献 [6] 所 采用 的 方法 只 
能 保证 载荷 的 释放 速度 限定 在 某 一 区 间 中 而 无 法 保证 释放 
速度 为 某 一 精确 值 , 且 未 对 弹簧 进行 空间 限制 ,因此 该 方法 
不 适用 于 需要 保证 精确 释放 速度 的 多 载荷 协同 编队 任务 。 

此 外 ,由 于 在 实际 分 离 过 程 中 不 可 避免 地 存在 阻力 和 其 
它 扰动 ,理论 分 析 无 法 准确 得 到 载荷 的 实际 释放 速度 ,为 保 
证 载荷 分 离 速 度 的 精确 性 ,还 需 对 分 离 过 程 进 行 Adams 运动 
学 仿真 计算 。 然 而 在 目标 速度 确定 的 情况 下 , 需要 试验 多 组 
数据 才能 得 到 满足 要 求 的 弹 段 性 能 参数 ,时 间 周 期 较 长 ,并 
且 设 计 结果 无 法 保证 质量 最 轻 的 设计 目标 。 因 此 为 减 小 航 
天 任务 中 的 运载 成 本 , 保证 载荷 的 顺利 编队 和 弹簧 的 力学 性 
能 ,对 多 载荷 分 离 释 放 弹 簧 进行 理论 分 析 与 仿真 计算 相 结合 
的 优化 设计 显得 尤为 重要 。 本 文 以 两 端 固定 、 不 磨 平 的 螺旋 
圆柱 弹簧 为 例 ,根据 性 能 约束 ` 几 何 约束 、 弹 自 可 靠 性 与 稳定 
性 的 要 求 ,建立 了 非 线性 约束 优化 模型 ,利用 Matlab 优化 工 
具 箱 求 得 弹簧 的 理论 最 优 解 , 将 理论 最 优 解 带 入 Adams 中 进 
行 运 动 学 仿真 ,得 到 理论 速度 与 仿真 速度 的 拟 合 关系 , 从 而 
对 于 航天 分 离 释放 任务 所 要 求 的 目标 速度 ,利用 拟 合 关系 得 
到 Matlab 中 所 对 应 的 理论 速度 ,进而 通过 Matlab 优化 工具 箱 
求 得 满足 特定 仿真 释放 速度 的 分 离 释 放 弹 簧 设计 参数 的 最 
优 解 。 


2 设计 原理 过 程 


2.1 设计 变量 

选取 三 个 独立 设计 参数 为 设计 变量 : 自 丝 直径 4d, 弹 段 中 
径 D, 弹 往 有 效 耸 数 n" , 故 令 设 计 向 量 X= [4d,D,n]” = 
[x ;x xs] 式 中 : D,d 为 连续 变量 ,n 仅 取 整 数 或 小 数 部 分 
为 0.5, 其 余 参 数 设 为 定 值 。 
2.2 目标 函数 

在 有 效 载荷 释放 任务 中 , 为 减轻 整体 质量 ,提高 运载 能 
力 , 以 弹簧 质量 M 最 小 为 优化 目标 。 可 得 


M = LpNDa Cn + 2)= i +2)N (1) 


4 4 
2.3 约束 条 件 
2.3.1 自由 长 度 


弹 段 自由 长 及 =pz+3d=(+3)d+56m 根 据 弹 段 安 
装 的 空间 要 求 得 : H, < 及, 即 
gi(X) =H-(x +3)x -x0 (2) 
式 中 : pi 为 自由 状态 下 弹簧 的 节 距 ,6 为 自由 状态 下 弹簧 自 
丝 之 间 的 距离 ,有 H 为 几何 空间 约束 下 所 人 允许 的 弹簧 最 大 
长 度 。 
2.3.2 装载 长 度 


与 弹簧 自由 长 类 似 ,弹簧 的 装载 长 度 刀 = (n+3)d + 5om， 
又 因 弹 簧 压 并 高 度 H, = (n+3) d, 可 得 
gl(X) =h -Cx + 0 (3) 
2.3.3 外 径 
因 弹 簧 外 径 不 得 超出 安装 空间 , 故 刀 +d<4, 其 中 4 为 
安装 空间 所 人 允许 得 最 大 外 径 。 可 得 
gi(X) =A-xl-x, 三 0 (4) 
2. 3.4 ”弹簧 指 数 


弹簧 指数 = 因 弹 簧 指数 减 小 会 导致 局 部 应 力 过 


大 。 另 外 ,弹簧 指数 较 大 或 较 小 时 ,加 工 性 也 成 为 问题 。 热 
轧 成 型 时 ,弹簧 指数 可 在 4 15 的 范围 内 选择 , 即 4<c<15， 
可 得 


2.3.5 纵横 比 
纵横 比 a = 分 ,纵横 比较 大 的 压缩 弹 和 在 轴 向 载荷 达 


到 一 定 程度 就 会 产生 侧 向 弯曲 而 失去 稳定 性 。 为 保证 使 用 
稳定 ,对 于 两 端 固定 的 压缩 弹 咎 ,纵横 比 一 般 选 为 : 0. 8 <a < 


H 
5.3 即 0.8< <5.3 即 


h— +3 
ge(X) =h+(6 es 上 四 0. 8x, 三 0(7) 
3 
h— 3 
g(X) =-h-(6, ee Die +5.3x, >0 
3 


(8) 

2.3.6 有效 省 数 
如 果 有 效 圈 数 为 3 以 下 , 弹 段 特 性 则 会 变 得 不 稳定 , 因 
此 应 将 弹簧 有 效 耸 数 设 定 为 3 以 上 , 即 ”>3, 结 合式 ( 3) 得 : 


4< 才 
即 
h 
gs (A) 0 (9) 
2. 3.7 应 力 约束 
弹 筑 中 的 应 力 
(10) 
式 中 : K 为 曲 度 系数 ,由 下 式 计算 = 二 1 + 全 ,为 工 
作 载荷 。 
最 大 工作 载荷 
Fe = 入 瓦 -有 (11) 


式 中 : 弹 咎 刚度 k= ,6 为 材料 的 剪 切 弹性 模 量 。 
则 最 大 应 力 


_ 8KDF。， _ CdCH, - “ 俐 -1 .0 | 


Dn 4c— 4 C 


s 

( 
式 中 : 7 为 弹 纂 材料 的 许 用 切 应 力 ,S, 为 许 用 安全 系数 , 因 
弹簧 要 求 精度 较 高 , 取 5, =1.3 ~1.7。 即 


B89 S, 
i ere ap ( 13) 
TTX2MX3 422 一 4 02 


2.3.8 性 能 约束 
为 保证 载荷 释放 时 达到 所 需 的 速度 ,根据 能 量 守恒 得 


上 -1 
7 (机 一 j02 = 


式 中 : mo 为 载荷 的 质量 ,m, 为 弹 咎 与 载荷 之 间 挡 板 的 质量 ， 
"为 载荷 的 释放 速度 ,w 为 摩擦 系数 ,7 为 载荷 所 受到 得 预 紧 
力 ,* 为 载荷 从 静止 到 释放 所 做 的 位 移 。 


(mo +m J) +uTs 


即 
tN 
gi0(X)= | Cx, + 3X + OX3 —h] 
160x2x3 
-Cm tm -an =0 ( 14) 


2.4 建立 数学 模型 
令 X= [dpD;a]” = [xxzx] ,得 到 以 圆柱 螺旋 压缩 
弹簧 的 质量 为 目标 画 数 的 非 线 性 约束 优化 模型 


minfCX) = mprix, Cx + DN 


1 
4 
(15) 


g(X) =0 (=1,2,.,9) 


gio( X) =0 
3 设计 计算 


3.1 设计 前 提 
将 载荷 在 地 面 实验 时 受到 的 重力 等 效 为 实际 分 离 时 的 
预 紧 力 , 故 预 紧 7 = mog。 在 载荷 分 离 过 程 中 ,不 考虑 除 摩擦 
阻力 之 外 其 余 阻力 对 载荷 运动 过 程 的 影响 。 优 化 计算 所 涉 
及 的 其 它 参 数 如 表 1 所 示 。 
表 1 弹 篱 设计 相关 参数 


自 参 数 参数 值 
材料 不 锈 弹簧 钢 
弹性 模 量 E/GPa 193 
剪 切 模 量 G/Gpa 71.6 
许 用 剪 切 应 力 7, /MP。 533 
密度 p/( kg/m’) 7.93 x103 
载荷 质量 mo /kg 10 


载荷 尺寸 /mm 100 x 200 x 300 


簧 装载 空 间 /mm 100 x200 


弹簧 参数 参数 值 
弹簧 与 载荷 之 间 的 挡 板 尺寸 /mm 3 x100 x100 
自由 状态 簧 丝 间 距 51 /mm 7 
弹簧 最 大 自由 长 H/mm 385 
挡 板 质量 mi /g 475.8 
滑 轨 摩擦 系数 | 0. 004 
重力 加 速度 g/( N/kg) 9.8 
安全 系数 5， 1.3 
容许 速度 误差 5% 
3.2 计算 过 程 
3.2.1 ”Matlab 优化 计算 


对 于 单 弹簧 设计 方案 ,为 保证 运算 速度 , 优化 设计 结果 
采用 Matlab 优化 工具 箱 中 的 有 效 集 算法 " , 对 于 不 同 的 分 
离 速 度 ,得 到 相应 的 设计 参数 ,如 表 2 所 示 。 

表 2 不 同 理论 分 离 速 度 对 应 的 设计 参数 


理论 分 离 速 ” 线 径 中 径 。 弹簧。 装载 力 弹簧 刚度 

度 vi(m/s) d(mm) D(mm) 着 数 n FF(N) k(N/mm) 
1 0.5 2.7 25,.1 13 103. 62 2.1598 
2 0.6 2.8 26.8 13 119.02 2.3136 
3 0.7 3.4 28.2 11 189. 42 4.7936 
4 0.8 3.6 30.2 11 213. 00 5.01 
5 0.9 4.3 36.8 10 259.11 6.1951 
6 1 4.3 33.9 10 314.65 7.72 
7 1.1 5 42.5 9.5 346. 25 7.978 
8 1.2 5 48.4 11 280.62 4.5161 
9 1.3 5.3 52.6 11 300. 63 4.5 
10 1.4 5.6 56.2 11 322. 81 4. 5535 
11 2 6 52 11 561.68 7. 13 


3.2.2 Adams 仿真 计算 
建立 如 图 1 所 示 的 设计 模型 。 


1 分 离 释 放 装 置 的 简化 结构 图 


将 表 1 中 的 装载 力 和 弹 纂 刚度 参数 带 入 Adams 模型 中 ， 
进行 仿真 计算 ,以 理论 分 离 速 度 v, = 0. 6m/s 为 例 ,得 到 如 图 
2、 图 3 所 示 的 载荷 位 移 曲 线 与 速度 曲线 。 

如 图 2 所 示 , 横 坐标 为 运动 时 间 t, 纵 坐标 为 载荷 质心 相 
对 于 分 离 装置 后 挡 板 的 位 置 s, 载 荷 初始 位 置 为 240. 5mm ,由 


Length/mm 


0.0 0.1 02 0.3 0.4 0.5 06 0.7 


Time/sec 


2 ， 理论 分 离 速度 v =0. 6m/s 时 ,载荷 的 时 间 一 位 置 曲线 


2 可 得 : 位 移 随 着 时 间 的 增 大 而 增 大 ,根据 分 离 装 置 的 设 
计 尺 寸 ,可 得 载荷 在 s =540. Smm 时 与 释放 装置 分 离 , 此 时 运 
动 时 间 t=0. 5195s。 
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3 “理论 分 离 速 度 v =0. 6m/s 时 ,载荷 的 时 间 一 速度 曲线 


如 图 3 所 示 , 横 坐标 为 载荷 的 运动 时 间 t, 纵 坐标 为 载荷 
的 速度 v, 根据 图 2 中 载荷 与 释放 装置 分 离 时 的 时 间 t = 
0. 5195s, 得 到 载荷 的 仿真 分 离 速 度 v, =0.533m/s。 以 此 类 
推 ,得 到 11 组 理论 分 离 速度 v, 及 其 所 对 应 的 仿真 分 离 速 度 
Vv; 如 表 3 所 示 。 
表 3 Matlab 理论 分 离 速度 与 Adams 仿真 分 离 速度 


vi/( mm/s) va/( mm/s) 相对 误差 
1 500 436 12. 84% 
2 600 533 11.17% 
3 700 613 12. 43% 
4 800 716 10. 51% 
3 900 816 9.31% 


Vi/( mm/s) v /1( mm/s) 相对 误差 
6 1000 906 9.45% 
| 1100 1000 9.09% 
8 1200 1127 6.08% 
9 1300 1229 5.48% 
10 1400 1330 5% 
11 2000 1912 4.41 


3.2.3” 拟 合计 算 

因 载 荷 在 释放 过 程 中 , 挡 板 与 载荷 之 间 存 在 接触 与 碰 
撞 , 分 离 过 程 存在 能 量 损失 ,所 以 仿真 释放 速度 小 于 理论 释 
放 速 度 。 由 表 3 可 得 ,利用 Matlab 计算 所 得 到 的 理论 值 与 
Adams 仿真 值 之 间 的 误差 在 速度 小 于 1. 4m/s 时 ,相对 误差 
均 大 于 5% ,超出 了 弹簧 设计 的 性 能 要 求 。 因 此 ,为 保证 弹 纂 
性 能 在 合理 的 误差 范围 内 ,将 表 3 中 的 数据 利用 Matlab 进行 
四 次 多 项 式 拟 合 ,得 到 理论 分 离 速 度 v 与 仿真 分 离 速 度 w 
之 间 的 拟 合 多 项 式 : 
oj = 0.1318v: - 0. 4469w + 0.4063w + 0. 9647w + 0. 0331 

(16) 

对 式 ( 16) 作 图 ,如 图 4 所 示 , 横 坐标 为 仿真 分 离 速 度 v, , 纵 坐 
标 为 理论 分 离 速 度 vi ,离散 点 均匀 分 布 在 曲线 两 侧 , 拟 合 精 
度 较 高 。 


vl/(m 
入 


In 一 


人 二 0 1 vi 


v2/(m/s) 


4 散 点 图 与 多 项 式 拟 合 曲线 


根据 式 ( 16) 中 拟 合计 算 的 结果 ,可 以 求 出 满足 特定 分 离 
速度 条 件 下 ,载荷 所 需 达 到 的 理论 释放 速度 ,进而 利用 Mat- 
lab 有 效 集 算法 反 求 出 弹簧 的 最 优 设计 参数 。 

3.3 ” 算 例 

在 多 载荷 协同 编队 探测 过 程 中 ,所 规定 的 载荷 释放 速度 
共有 三 组 ,分 别 为 :0. 593m/s、0.799m/s、0.92m/s。 原 始 的 设 
计 参 数 如 表 4 所 示 。 


表 4 分 释放 弹簧 的 原始 设计 参数 


组 目标 速度 线 径 中 径 ”弹簧 质 ”真实 仿真 速度 
vo/(mm/s) dmm D/mm 量 m/g va3/(mm/s) 
1 593 5.01 63 390 593. 6 
2 799 5.59 61.86 430 799.3 
3 920 5.89 60.44 450 919.6 


以 目标 速度 v。= 0.799m/s 为 例 ,假定 有 效 载荷 最 终 的 
释放 速度 与 Adams 的 仿真 结果 相同 , 即 另 v, = vw, 根据 式 
( 16) 得 载荷 理论 分 离 速 度 v, = 0. 889m/s, 利 用 Matlab 优化 
工具 箱 中 的 有 效 集 算法 得 到 : 装载 力 F = 237. 1825N, 刚 度 k 
=5.2915N/mm,; 线 径 d =3.8mm, 中 径 D =31.9mm, 弹 得 有 
效 耸 数 n = 11. 将 上 述 参 数 带 入 Adams 模型 中 进行 仿真 ,得 
到 载荷 的 实际 仿真 分 离 速度 v, =0. 805m/s. 同 理 对 于 其 它 两 
组 目标 速度 ,得 到 如 表 5 所 示 的 速度 和 质量 对 照 表 。 

表 5 载荷 分 离 速度 与 质量 对 照 表 


目标 速 ” 理 论 速 仿真 速 相对 优化 原始 质量 缩 

度 w。 度 v 度 v 误差 质量 质量 减 率 
1 593 671.1 606.9 2.3% 73.4 390 81.2% 
2 799 889.0 805.0 0.8% 117.3 430 72.7% 
3 920 1010.9 927.4 0.8% 143.5 450 68.1% 


由 表 5 可 得 ,实际 仿真 速度 与 规定 速度 的 误差 均 在 5% 
以 内 ,满足 弹 咎 设计 的 性 能 要 求 , 且 弹 簧 质量 有 明显 的 缩减 ， 
满足 优化 设计 的 要 求 。 


航天 任务 的 特殊 性 要 求 弹 簧 采用 最 优化 设计 ,并 且 载 荷 
在 释放 过 程 中 ,不 可 避免 地 存在 能 量 耗 散 , 理论 分 析 无 法 准 
确 得 到 载荷 的 实际 释放 速度 。 同 时 仿真 计算 存在 收敛 速度 
慢 ,时 间 周 期 长 ,无 法 进行 优化 设计 等 缺点 。 本 文 利用 最 优 
化 理论 ,建立 了 关于 弹簧 推进 释放 速度 的 理论 分 析 模 型 , 然 
后 利用 Matlab 拟 合计 算 与 Adams 仿真 计算 ,建立 理论 释放 速 
度 与 仿真 释放 速度 的 关系 式 , 再 利用 Matlab 优化 工具 箱 , 反 
求 出 满足 设计 要 求 的 弹簧 最 优化 设计 方案 。 本 文 的 设计 过 
程 具有 运算 速度 快 ,精确 度 高 ,参数 优化 等 优点 ,减少 了 仿真 
次 数 ,提高 了 释放 速度 的 精确 性 , 在 以 弹簧 作为 推进 装置 的 
其 它 机 构 中 也 具有 重要 的 参考 价值 。 
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